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t^ Abstract 

i— —i This paper discusses a constrained gravitational three-body problem with two of 

the point masses separated by a massless inflexible rod to form a dumbbell. The non- 
►>. integrability of this system is proven using differential Galois theory. 

On 

^ l Equations of motion, symmetries and reduction 

^O 

_t Considered is the gravitational three-body problem with a single constraint. Three point 

masses, m\, m 2 and m 3 move in a plane under mutual gravitational interaction. Masses 

m 2 and m 3 are connected by a massless inflexible rod of length / > to form a dumbbell. 

L* A pictorial description of the problem is given in Figure IT] 

• ^h The Lagrangian of the system has the form 

r> 111 

L = -mi||i-i|| 2 + -m 2 \\i 2 \\ 2 + -m 3 ||r 3 || 2 - U(x lr X2,x 3 ) r (1.1) 

where Xi, r 2 and r 3 are the position vectors of the respective masses, and 

Gm 1 m 2 Gm\m 3 Gm 2 m 3 

U(r 1 ,r 2 ,r 3 ) = - T n — Ti n" ~ Ti if' ( 1 - 2 ) 

II r l _ r 2|| || ti — 1*3 1| || r 2 — 1*3 1| 

is the potential, r d is the radius vector of the centre of mass of the dumbbell, and 

r 3 ~~ r2 n n\ 

e := ti rr/ (1-3) 

1 1 r 3 — r 2 1 1 

is a unit vector along the dumbbell. It follows that 

r2 = r d + fi 3 le, i 2 = id + fait, (1-4) 

r 3 = r d - fi 2 le, x 3 = r d - fi 2 le, (1.5) 




Figure 1: A dumbbell and a point mass in the plane. Point O denotes the origin of an 
inertial frame. Point C is the system centre of masses. Axis OX\ of the moving frame is 
always parallel to the dumbbell. 



where 



n- 



m d 



n v - 



m 3 

m d 



m d '■= m 2 + OT3. 



(1.6) 



The system has five degrees of freedom. The configuration space of the point m\ is M. 2 , 
and the configuration space of the dumbbell is M 2 x S . A configuration of the system is 
fully specified by ti, r d , and e. In these coordinates, the Lagrangian is as follows. 



where 



I:=l 



2 
m d 



mi ri 



1 1 ~ 

I 2 + 2 m d||rd|| 2 + 2 J II^II 2 ~ U(ti,t d ,e), 



U(r lf r d/ e) = U{t\, r d + }i 3 le, r d - ]i 2 le) 



(1.7) 



(1.8) 



The system possesses natural symmetries that are exploited to reduce the dimension of 
the system. Amongst these is translational symmetry. This is manifest by defining r : = 
r d — 1\ as the vector between the mass m\ and the dumbbell. A configuration of the system 
can be thus described by r, e, and r s defined by 



As 



m s r s = mm + m d r d/ 
m d 



nu 



t, 



nu 



*d 



m\ + ni2 + WI3. 
nil 



nu 



t, 



the Lagrangian ( |1.7| is written in the following form. 



-m r \\r\ 



-I\\e\ 



-m s ||r s || 2 - W(t,e), 



(1.9) 
(1.10) 

(1.11) 



where 



is the reduced mass, and 



m r 






W(r,e) 



Gniini2 
r + p 3 le\ 



Gm\m 3 
r — jiile\ 



(1.12) 



(1.13) 



The components of r s are cyclic coordinates, and the motion of the centre of mass separates 
completely. Therefore, the term m s ||f s || 2 /2 is removed from the Lagrangian ( 1.11) . 

It is convenient to introduce dimensionless variables, taking / as the unit of length, m T 
as the unit of mass, and 

T:=t/— , (1.14) 




as the unit of time. Setting / = m r = T = 1, the Lagrangian ( 1.11 1 reduces to 

L = ^l|r|| 2 + ^||e|| 2 -W(r,e), 
where 



W(r,e) 



¥■-- 



m 3 



(1.15) 



(1.16) 



|j (1 + f/)r + J we|| || (1 + }i)t — e||' ' m 2 

The reduced system still has a symmetry: it is invariant with respect to the natural action 
of the group SO (2, R). In fact, the dynamics are oblivious to the orientation of the inertial 
frame. This symmetry can be used to reduce the dimension of the configuration space 
by one. This is achieved by describing the dynamics in a rotating frame in which the 
dumbbell is at rest. The transformation from the inertial frame to the rotating frame is 
given by 

r = AR, AeSO(2,R). (1.17) 

An additional assumption, that e is parallel to x-axis of the rotating frame, allows for a 
coordinate representation of the transformation as follows. 



Then 
and 



cos cp — sin cp 
sin cp cos cp 

e = AE, 



and e = [cos <p,sin<p] T . 

E=[1,0] T , 



<p[-sin<p,cos<p] T , A = <pAJ, J 



-1 

1 



The Lagrangian expressed in the coordinates R = [Xi, X2J T , and cp has the form 

L = ^||R|| 2 + <pR T JR + ^<p 2 ||R|| 2 + \l<P 2 ~ V(R), 



where 



V(R) := 



/' 



(l + }i)R + }iE\\ ||(l + p)R-E| 



(1.18) 

(1.19) 
(1.20) 

(1.21) 
(1.22) 



The generalised momenta P := [Pi, P2] and P v , are given by 

P:= |i = R + <pJR, P<p := ^ = <p(l|R|| 2 + I) +R T JR- (1-23) 

Coordinate q> is cyclic, so P,p is a first integral of the system. Thus, the Hamiltonian of the 
system is written as 

H = l\\n 2 + \oc (7 - P T JR) 2 + V(R), (1.24) 

where a. = II I, and 7 is a fixed value of the first integral P«, corresponding to the cyclic 
variable cp. Hence, the equations of motion have the form 

^R =P - a (7 - P T JR) JR, 

^P = - a ( 7 -P T JR)jP-(l + F)# ± 4 R ± i ^- (125) 

df V i )i \ p; || (1 + ^)R + ^E|| 3 Kl - A0 > 

. (l + u)R-E 



(1 + ^)R-E|| 3 ' 



2 Problem 

Numerical experiments suggest that the considered system is not integrable. The complex 
behaviour of the system is apparent from the Poincare cross-section. For fixed values 
of the parameters and the energy the equations ( |1.25[ ) are integrated numerically by the 
Burlish-Stoer method 0. The cross-section plane is specified as Xi = 0, and the cross 
direction is chosen Xi > 0. The coordinates of the cross-section are (X2, Pi). The cross- 
sections presented in the subsequent figures are obtained for parameters y. = 1 and oc = 
i_ 
ex- 

Figure [2] is a cross-section corresponding to energy e = —0.1 and 7 = — 1. Chaotic 

behaviour is evident. Two cross-sections for 7 = are shown in Figure [3] and Figure |4l 
Both of them attest to the non-integrability of the system. The problem is proving the 
non-integrability of the system analytically. It is attractive to try using the differential 
Galois approach to the problem of integrability (see e.g. |4J). However, as is explained in 
|U, application of this theory directly to the system ( 1.25 1 is invalid as the Hamiltonian 
function ( |1.24| | is not single-valued. 

Proposed in this paper is a solution addressing this deficiency, based on an extension 
of the system into a larger phase space. Coordinates of these additional dimensions are 
denoted by Ri and Ri- This approach further requires a definition of the time derivative 
consistent with system ( 1.25) . Following this, coordinates R\ and R2 are considered as 



lengths of vectors 

Ri:=(l + f/)R + ^E, R 2 :=(l + p)R-E. (2.1) 

That is, equations G\ = G2 = 0, where polynomials G\, G2 £ C [Xi, X2, R\, R2] given by 

G 1 :=Rf-[(l + /.)X + M )] 2 -(l + /.) 2 Y 2 = R?-R[R 1 , 1 
G2:=R 2 -[(l + ^)X-l)] 2 -(l + /Y) 2 y 2 = Ri-R 2 r R 2 . J 




Figure 2: Cross-section for energy e = —0.1 and 7 = — 1 

defines R\ and R2 as algebraic functions of (Xi, X2). 

The extended phase space is C 6 with coordinates Z := [Xi, X2, Pi, P2, Ri, R2] • The 
Hamiltonian ( 1.24[ > expressed in these coordinates reads 



K = K(Z) := hn 2 + 1 ct( 1 - P T JR) 2 - — - £. 



As 



d7 R < R 



R ^(1 + ,) R ; R 



K, ' 



* = 1,2, 



(2.3) 



(2.4) 



system ( 1.25| l can be extended to the following one. 
d 



df 



Z = U(Z), U(Z):=J(Z)K'(Z), J(Z) 



2 E 2 


2 ] 


E 2 2 


-A 


2 A T 


2 J 



where 



A=(l + p) 



Ri R2 

R t 'R 2 



K'(Z):=(V Z K(Z)) T . 



(2.5) 



(2.6) 



It can be shown that matrix J(Z) defines a Poisson structure on C 6 . The corresponding 
Poisson bracket is given by the following formula 



{F,G}(Z)=F'(Z) T ](Z)G'(Z), 



(2.7) 



1 - 



CM A 

Oh U 






#"V' 



•""'^7/- ■> -y? 



■>-■-■ 






-4 



; : i ; '((g)) { (°))i & -> Ml ^SP 



^■t 



5^>^S: 











x 2 



Figure 3: Cross-section for e = —0.2 and 7 = 

where F and G are two smooth functions. The Poisson structure given by J (Z ) is degener- 
ated and has rank 4. It is demonstrable that polynomials G\ and G 2 are Casimir functions 
of this structure, whose common levels are symplectic manifolds. In particular, 

MJ:={ZgC 6 | Gi(Z) = G 2 (Z) = 0} / u G (0,1/2]. (2.8) 



is a symplectic manifold. System ( |2.5| , when restricted to Mi is equivalent to ( |1.25[ ). 

In the setting thus described the main result of this paper can be formulated in the 
following theorem. 

Theorem 2.1. For 7 = and ]i G (0, 1/2], system ( |2.5| l, when restricted to the symplectic leaf 
Mt, is not integrable in the Liouville sense with meromorphic first integrals. 



3 Particular solution and variational equation 

If 7 = 0, then equations ( |2.5[ ) have an invariant algebraic manifold, 

X := {Z G C 6 | X 2 = P 2 = 0, Gi(Z) = G 2 (Z) = } . 
On this manifold it can be taken that 

R 1 = (l + 1t)X 1 + n, R z ={l+jt)X 1 -l, 
so the system ( |2.5| l, restricted to N, reduces to 

P = -(p + l)[ 2 , + 1 = 



X = P, 



(3.1) 
(3.2) 

(3.3) 



1 - 



Oh U r 



fm 



/l -~f V v f \ 





x 2 



Figure 4: Cross-section for e = —0.15 for system with 7 = 

where (X,P) := (Xi,P\) £ C 2 . This is a Hamiltonian system with one degree of freedom. 
Thus the manifold N is foliated by phase curves r e on energy levels Km = e, that is 



r t cM,:= (X,P)GC 2 e 



/' 



(l + p)X-l {l + }i)X + ii 



(3.4) 



For a generic value of e the level M e contains three phase curve. If e is real, two or three 
of these levels have a non-empty intersection with the real part of the phase space. In 
further consideration it is assumed that e > 0, and the chosen phase curve T e contains the 
half-line (1/ (1 + yi), 00) x {0} C C 2 . In the phase space C 6 this curve is given by 



r e :={zex| (Xi,l*L)er e } 



(3.5) 



Application of the Morales-Ramis theory necessitates the linearisation of equations ( 2.5 1 
along T e . It has the form 



dt 



Vz, V= f^( Z )' Zer <- 



(3.6) 



This system has three first integrals 

fc(z):=VK(Z)-z, gi(z):=VGi(Z)-z / g 2 (z) := VG 2 (Z) • z, ZeT, (3.7) 
On the level gi(z) = £2(2) = 0, the following equalities are satisfied 

R 1 r 1 = (l + ii)[(l + ii)X 1 + fi\x 1 , R 2 r z = (l + ii)[(l + ii)X 1 -i\x 1 , (3.8) 



where it is assumed that z = [x\, x 2 , p\, p 2 ,r\,r 2 ] ■ On T e , coordinates (R^Rz) can be 
expressed by X\, 

Ri = (1 + /<)X! + ji, R 2 = (1 + p)X! " 1- (3-9) 

Hence, on the level gi (z) = g 2 (z) = 0, variables (j\, r 2 ) can be eliminated from the varia- 
tional equations (|3.6|. The reduced system of variational equations has the form 



(3.10) 





Ml 




"0 1 





0" 
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V\ 
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V\ 


dt 


u 2 







-aXP 


1 + aX 2 




u 2 




Pi. 




.0 
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ocXP_ 




y 2 



where 



A=2{ } i + l) 



¥ 



+ 



[(i + ¥ )x-iy [(i+rix+tf 



(3.11) 



B 



-ocP 1 



A, 



and {X,P) S T e . This system splits into two subsystems. Variables {u\,V\) describe varia- 
tions in the invariant plane N. In fact, [u\, V\) is a vector tangent to N at point (X, P), hence 
the subsystem in variables {u\, V\ ) is called tangential. The second subsystem, correspond- 
ing to variables (u 2r v 2 ), yield the normal variational equations. The latter subsystem can be 
expressed by a single second-order equation, 



where u = u 2 , and 



a :- 





u + aii + bu = 0, 


locXP 


IccP 1 


1 + ocX 2 ' 


1 + ccX 2 



(3.12) 



fl + aX 2 ) A 



Changing the independent variable by 

t — >z = {\ + }i)X 



(3.13) 



transforms equation ( 3.12} into an equation with rational coefficients. This follows from 



d_ 
dt 



dt 2 



(i) 



,d 2 



+ z- 



dz' dt 1 w dz 2 dz' 
z = (1 + ji)x = (1 + ¥ ) Pl , (z) 2 = (1 + u) 2 V \ = 2(1 + ¥ ) 2 [ e 

Y 1 



z + u 



Z= (l + ji)p 1 = -(fi + 1 
where the energy integral 



+ 



-X 1 



(z-l)2 (z + ¥ ) 



¥ 



(1 + p) X - 1 (1 + p) X + n 



has been used. 



Resultantly, equation ( 3.12| | transforms into 



X" + p(z)X' + q(z)X = 0," = — , X = X lr 



(3.14) 



with 



p = - [(-1 + z) 2 (l + z 2 (5 + 4ez)cc) + (2 + z(-4 + 3z + (4 + z(-8 + 8e(-l + z) 2 
+5z 2 ))tx))p + (1 + 3z 2 + 2z(2e(-l + z) 2 + z(-4 + 5z))oc)p 2 + (1 + z(4 + z - 4a 

+5za))p 3 + 2(1 + z)p 4 + p 5 ] / [2(z - l)(z + p)((z - 1)(1 + ez) + (e(z - 1) + z)y. 

V)(^ + (i + f) 2 )], 

<? = [4(z - l) 2 a(l + p)(z + }i) 2 ((z - 1)(1 + ez) + (e(z - 1) + z)p + p 2 ) + ((z - l) 3 
+ (1 + 3(z - \)z)y. + (3z - l)p 2 + p 3 )(z 2 a + (1 + p) 2 ) 2 ] / [2(z - 1) 2 (1 + jt) 
•(z + p) 2 ((z - 1)(1 + ez) + (e{z - 1) + z)p + p 2 )(z 2 a + (1 + p) 2 )] 

Making the next change of variables, 



X = w exp 
( 3.14| | is transformed to the canonical form, 



p(s)ds 



zo 



1 1 

w" = r(z)w, r(z) = -q(z) + -p'{z) + |p(z) 2 , 



(3.15) 



(3.16) 



with coefficients given as 
P 



M' 



2l2 



M = 16(z - \ff{z 2 + /6 2 ) 2 (z + ^) 2 [(z - 1)(1 + ez) + (e(-l + z) + z)jt + p 2 } 2 , (3 17) 



li 



p=j:p,z u -. 



i=0 



Here /3 2 = (1 + }i) 2 / a and coefficients p, of P are given in the Appendix in ( A.1) . This 
equation has seven singularities, 

zi =0, z 2 = -p, z 3/4 = ±ij8, 

e - 1 - (1 + e)p ± V(l + j/)((l + e)2 + ( e - 1)2^) 



25,6 



2c 



and Z7 = 00, providing energy is chosen such that 

r (l±i J S)(±i + £Tip)(£±ip)(l + p) (v^ii) 2 



e^ 



(1 + ^2)^2 + ^) 

p + p 3 ±i(i + p 3 + p 2 (i + p))\p\\ 



p + l 



(3.18) 



(1 + P)<JP + 1 *) 



Points Zi, . . . ,Ze are regular poles of second order and the order of infinity is l provided 
that E(l + p) 7^ 0. Here, the order of a singular point is defined as in [2j. 



The differences of exponents at the singularities in C are as follows. 

A 2 = 



P 2 ' 




A3 



2, A 5 



1 
2' 



(3.19) 



It is evident that the differences of the exponents at singularities Z3 and Z4 are integer, 
and thus, in local solutions around these points, logarithmic terms may appear. Appli- 
cation of the method described in Appendix B verifies the presence of such terms. For 
singularity Z3, solutions of the indicial equation are oc\ = 3/2 and a.2 = —1/2. The rele- 
vant one is a = ct\ = 3/2. The expansion of r{z — Z3) 2 , according to (B.4 1, gives coefficients 
ro, t\, and Y2- Then/1,/2 and g\,gi are determined from ( |B.6| and ( B.12| respectively. Since 
s = oc\ — 0.2 = 2, the coefficient g2, which multiplies the logarithm, must be found, and it 
is given as 



c?2 



i(l + f /)((i + J 6) 2 + ^-2i^-^ 2 ) 



8^i + ^-iy)(l + e^ + (e-y)y-i^l + e(-l + y) + y)Y 
Examination of the real and imaginary part of g2 yields the following conditions. 

/3(1 + ji)[l + yS 4 (l + y) + 2£ 2 (-l + }i)(-l - ley + y 1 ) 
+ y{-\ - 2e(-l + y)(l + y 2 ) + y(-l + y(-l + (-1 + y)y)))] = 0, 
(l + y)[y(l + 3p 2 + (-l + y)y)(-l + y 2 )+e(f5 6 + y 2 (-l + y-y 2 ) 
+ p i (2 + y(-3 + 2y))+p 2 (l + y(-l + y(-l + (-l + y)y))))\ 

This system has two solutions satisfying y > 0, ft 2 > and e£R. 

1 + y 



(3.20) 



0. 



^ 



¥ 



y 2 , e 



V 



and j6 = y — 1. 



(3.21) 



The condition precluding a logarithm in a local solution around Z4 is g| = 0/ where * 
denotes complex conjugation. The corresponding solutions are the same as those given 
previously. 

In the first solution in ( 3.21) , condition 



-1 + y-y 2 



(3.22) 



gives 



mo 
nil 



y{y + l) 



<0, 



y 2 + l 

■ 0. The second solution in ( 3.21) also yields 



with the only non-negative solution m.\ 
mi = only 

If conditions ( 3.21) are not satisfied, the two linearly-independent local solutions W\ 
and W2 of ( |3.16) in a neighbourhood of z* = Z3 or z* = Z4 have the following forms. 

z,) + {z-z,f- 2 h{z) l 



w\{z) = (z — z t .) a f(z), W2(z) = W\(z) ln(z 



(3.23) 



where f(z) and h(z) are holomorphic at z* and /(z*) 7^ 0. The local monodromy matrix, 
which corresponds to the continuation of the matrix of fundamental solutions along a 
small loop encircling z* counter-clockwise, gives rise to a triangular monodromy matrix, 







-2m 
-1 



10 



(For details, see 0.) A subgroup of SL(2, C) generated by a triangular non-diagonalizable 
matrix is not finite, and thus the differential Galois group is not finite either. Moreover, 
the differential Galois group Q of this equation is not any subgroup of the dihedral group 
because such subgroups contain only diagonalizable matrices. Thus, Q is either the full 
triangular group or SL (2, C ) . But since the order of infinity is one, the necessary condition 
that Q is the full triangular group is not satisfied (see Lemma [B.2| . This means the that 
the differential Galois group of ( |3.16 l is SL(2, C), with a non-Abelian identity component 



equal to the whole group SL(2, C). 

Regarding exceptional values of energies given in ( 3.18| |: the condition that the first 



one is real gives /3(1 + }i 3 + fi 2 (l + }i)) =0, which further suggests the equality ( |3.22| l. Its 
only physical solution is mi. The second expression in ( 3.18| l is never real and the third 
one is real when ( 3.22| | holds only. 
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(A.1) 



A Appendix 

Explicit form of coefficients p, in r given by §3.17) are as follows. 

po = -8e(l + //), vi = "8(1 + p)(l +4e(-l + ?0 + ?<)]/ 

p 2 = -8(l + //)(-4 + 4// 2 + e(6 + J 6 2 + 2//(-5 + 3//))) / 

p 3 = -48e 2 i 6 4 - 16e(-l + // 2 ) (2 + 20 2 + p(-3 + 2//)) - (1 + /^) 2 (11 J 6 2 

+ 24(2 + //(-3 + 2//))), 

p 4 = -4(48e 2 J 6 4 (-l + y) + (-1 + p)(l + ^) 2 (8 + ll^ 2 - 6// + 8// 2 ) 

+ 2e(l + //)(1 + 120 4 + y6 2 (6 + //(-7 + 6//)) + ?*(-5 + ^(5 + (-5 + //»))), 
p 5 = -2(-4(l + (e- p)n) (-1 + // 4 ) + 3j8*(9(l + }if + 64e(-l + // 2 ) 

+ 16e 2 (3 + //(-8 + 3//))) + i6 2 (16e(-l + // 4 ) + 3(1 + //) 2 (11 + p(-13 + 11/*)))), 
p 6 = -4 J 6 2 (48e 2 J 6 2 (^ - 1)(1 + (// - 5)/*) + (// - 1)(1 + //) 2 (11 + 540 2 

+ //(3 + 11//)) + 2e(l + ji)(1 + 3^(24 + //(-59 + 24//)) 

+ f*(-l + M-5 + (f-l)f*))))/ 

p 7 = -0 2 (11 + 30 4 (1 + //) 2 + 12 J 6 2 (4e(-8 + 33// - 33// 3 + 8// 4 ) + (1 + //) 2 (27 

+ //(-59 + 27//)) + 4e 2 (l + (-2 + //)//(8 + (-14 + //)//))) + //(32 + //(-42// 

+ // 3 (32 + 11//) - 32e(-l + // 2 )))), 

p 8 = -4 i 6 2 (3£ 4 (-l + //)(l + //) 2 - 2(1 + (g -//)//)(// + // 4 ) 

+ 2 J 6 2 (-24e 2 (-l + //)//(l + (-5 + //)//) + 27(-l + //)(1 + //) 2 (1 + (-3 + //)//) 

+ 2e(l + //)(6 + //(-76 + //(157 - 76// + 6// 2 ))))), 

p 9 = -2j6 4 (27 + 3/3 2 (l + //) 2 (3 + //(-5 + 3//)) + //(-176 - 140e(-l + 4// - 4// 3 

+ // 4 ) + 48e 2 //(3 + //(-8 + 3//)) + //(4 + //(414 + //(4 + //(-176 + 27//)))))), 

pio = -4 J 6 4 (3 J 6 2 (-1 + // 2 - // 3 + // 5 ) - 2//(l + (e - //)//)(-H + //(24 

+ 24e(//-l) + (24-ll//)//))), 

pu = -3/3 4 (16// 2 (l + (e - //)//) 2 + £ 2 (1 + // 3 ) 2 ) 

B Second-order differential equation in reduced form with rational coeffi- 
cients AND ITS DIFFERENTIAL GALOIS GROUP 

Consider the second-order linear differential equation in reduced form with rational co- 
efficient 

w" = r(z)w, r(z) G C(z). (B.l) 

A point z = c G C is a singular point of this equation if it is a pole of r(z). A singular 
point is a regular point if at this point (z — c) 2 r(z) is holomorphic. 

Assume that c is a regular point and look for local solutions of this equation in a neigh- 
bourhood of this point. For simplicity, assume that c = 0. Look for a solution of the 
following form (see e.g. [6]), 

oo 

w = zj(z), /(z) = l + X>*', (B.2) 

12 



where oc,fi for i S N are constants to be determined. Substituting in (|B.1[) yields 



,«-2 



r{z)z a \l + YLfi zi ) =°- ( R3 ) 



«(« - 1) + £]0 + z)0 + i ~ l )fi z 

1=1 

Multiplying by z 2 , apply Taylor's formula for the analytic function z 2 r{z) 

CO 

z 2 r(z) =r + J] nz\ 



(B.4) 



i=\ 



and equating to zero the coefficients of successive powers of t, obtain the sequence of 



equations determining unknown coefficients in series (B.2|. The first of these equations, 
called the indicial equation and determining a, has the following form. 



a(« — 1) — ro = 0. 
From the remaining equations / is evaluated. 



h 



2a' 



'i-l 



fi 



i(i + 2«-l) \J>*-i + ri )' fOT ^ Z 



(B.5) 



(B.6) 



Let ct\ and oti be the roots of the indical equation. If the difference s = oc\ — ot^ is not an 
integer, then equation ( |B.1| | has two local solutions of the form 



w 1 = z*i 1 + ^0 , w 2 = z^[l + £//z ! 



(B.7) 



i=\ 



i=\ 



In the case where s is a positive integer or zero, the series Wi may not exist or coincide 
with w\, and a second independent solution must be constructed as 



w 2 {z) = w\(z) 



dx 



wi(xy 



wi(z) 



\-2cg ^ 



/(*) 



wi(z) I z 6 ' 



f( x y 



(B.8) 



In the last equality the relation 2«i = s + 1 is applied. 
Denoting 

-| oo 

the second solution is written as 

w 2 {z) = w x (z)g s In z + z tt2 



+ E^ )> 



(B.9) 



(B.10) 



j=l 



where /z, are constants. 

Coefficients g,- are given by 



Si = -2/i, 



L ?2 



-2/ 2 -2/^! -A 2 = -2/2 + 3/2, 



and 



i-l t-2 J-/-1 

gi = -2fi-E( 2 $H+fH)fj-Lfi E /*#-/-*' 

;=1 ;=1 fc=l 



(B.ll) 
(B.12) 
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for i > 3, while coefficients hi are defined as follows. 

r h gi 



lh 



s s — 1 

s s 



Si 






A 

s 

fs+1 



;=i 



s-l 



; 



fs+i-j 



fi , gi 
s i — s 



7=1 
1 gj - [ 



s J ;=s+i / s 



fH 



for z = 1, 

for 1 < i < s, 

for z = s, 

for z = s + 1, 

for z > s + 2. 



(B.13) 



7=1 



;=s+l 



The differential Galois group of equation ( |B.1[ ) is a subgroup of SL(2, C). Classification 
of these subgroups is well known. The following lemma shows that the differential Galois 
group determines the form of solutions (see e.g. ©SI). 



Lemma B.l. Let Q be the differential Galois group of equation ( 3.16 1. TTzez-z one of four cases may 
occur. 

Case I G is conjugate to a subgroup of the triangular group, 



r 



a b 
a- 1 



a E C*,b G C}text; 



in this case equation ( |3.16| l has an exponential solution of the form y = Pexp J to, where 
P G C[z]andco G C(z). 



Case II Q is conjugate to a subgroup of 



V" 



c 
c- 1 



c G C* \ U 



cGC* text; 



in this case equation ( 3.16| l has a solution of the form y = exp J to, where co is algebraic 
over C(z) of degree 2, 



Case III Q is primitive and finite; in this case all solutions of equation ( |3.16| l are algebraic, thus 
y = exp J co, where co belongs to an algebraic extension ofC(z) of degree n = 4, 6 or 12. 

Case IV Q = SL(2, C) and equation ( |3.16| l has no Liouvillian solution. 



Kovacic |21 formulates the necessary conditions for the respective cases from Lemma B.l 
to hold. 

Write r(z) G C(z) in the form 

s(z) 



r(z) 



t(z)' 



s(z), t(z) G C[z], 



w^here s(z) and t (z) are relatively prime polynomials and £(z) is monic. The roots of £(z) 
are poles of r(z). Denote E' := {c G C | f(c) = 0} and E := 1/ U {°°}. The order ord(c) of 
c G E' is equal to the multiplicity of c as a root of t (z), the order of infinity is defined by 

ord(oo) := degf — degs. 
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Lemma B.2. The necessary conditions for the respective cases in Lemma^ETT are the following. 



Case I. Every pole ofr must have even order or else have order 1. The order ofr at oo must be even 
or else be greater than 2. 

Case II. r must have at least one pole that either has odd order greater than 2 or else has order 2. 

Case III. The order of a pole ofr cannot exceed 2 and the order ofr at oo must be at least 2. If the 
partial-fraction expansion ofr is 

then Aj = \/l + 4a ; - e Q,for each i, ^ bj = and if 
then ^1+4^ G Q. 
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